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Abstract 

Exciton-polariton condensates are considered as a deterministic source of bright, coherent non- 
Gaussian light. Exciton-polariton condensates emit coherent light via the photoluminescence 
through the microcavity mirrors due to the spontaneous formation of coherence. Unlike con- 
ventional lasers which emit coherent Gaussian light, polaritons possess a natural nonlinearity due 
to the interaction of the excitonic component. This produces light with a negative component to 
the Wigner function at steady-state operation when the phase is stabilized via a suitable method 
such as injection locking. In contrast to many other proposals for sources of non-Gaussian light, 
in our case, the light typically has an average photon number exceeding one and emerges as a 
continuous wave. Such a source may have uses in continuous-variable quantum information and 
communication. 

PACS numbers: 71.36.+c,42.50.Gy,42.50.Ex 
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Quantum optics represented by continuous variables [l-^ has been shown to be a useful 
approach to quantum information processing, complementary to the traditional discrete- 
variable methods using qubits. While simple, elementary quantum protocols such as tele- 
portation [6] can be realized using Gaussian continuous-variable (CV) states and Gaussian 
operations, any more advanced application such as universal quantum computation would 



require some non- Gaussian elements 
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-9|. Moreover, even non-Gaussian states may be in- 
sufficient as resources for quantum computation: trivially, those which are expressible as 
mixtures of Gaussian states; however, as well those which are not represent able as mixtures 



of Gaussian states and yet possess a positive Wigner function [10|, (for the discrete ana- 
logue, see [I4]). As a consequence, the occurrence of a negative Wigner function can be seen 
as a prerequisite for a potential quantum mechanical speed-up [13] (for related results on 



the discrete case, see 



14j). While such iiegativities are necessary, even only small instances 



of it are sufficient for universality 
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Unfortunately, since there are no media 



with strong enough nonlinearities, current experimental approaches to creating states with 
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18| . only conditionally 



preparmg 



a negative Wigner function are highly probabilistic 
such quantum states (those few deterministic schemes using solid-state emitters [19| would 
still not be capable of producing bright states). A more desirable way of generating non- 
Gaussian light is a more "plug-and-play" style device, where the device is simply switched 
on and bright non-Gaussian light emerges deterministically. 

In this paper we describe a deterministic method of producing non-Gaussian quantum 
states of light with a negative Wigner function using exciton-polariton condensates. Exciton- 
polaritons are bosonic quasi-particles in semiconductor microcavities, corresponding to a 
coherent superposition of an exciton and a cavity photon [20] . Exciton-polaritons have been 



observed to undergo a dynamical form of Bose-Einstein condensation |21l . |22| . One of the 
attractive features of this system is that the state of the polaritons inside the semiconductor 
can be directly imaged using the leakage of photons through the imperfect microcavity 
mirrors. Thus the spontaneous buildup of coherence, one of the signals of Bose-Einstein 
condensation, can be directly measured simply by analyzing the properties of the emerging 
light. 

One of the principle differences between an exciton-polariton condensate and a laser is 
the presence, or lack respectively, of strong coupling between the excitons and photons. 
Strong coupling leads to polariton-polariton interactions, giving much richer physics than 
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its non-interacting counterpart. These interaction effects also lead to phenomena such as su- 



perfluidity of the polaritons 



23|. In this paper, we take advantage of the polariton-polariton 



interactions to propose a device for producing non-Gaussian light. The key technological 
advance that has allowed for this is the availability of high-Q cavities where polariton life- 
times of the order of ~ lOOps are now possible 2J]. In addition, exciton-polaritons may be 
confined spatially using trapping techniques such as metal deposition on the surface of the 
microcavity [25] , increasing the self-interaction energy U. This gives the possibility of energy 
scales of the self-interaction U and the cavity decay ^70 to be of the same order, leading to 
the creation of non-Gaussian light in the steady state. One of the attractive features of this 
method is that the light is produced completely deterministically and continuously while 
it emerges from the microcavity. This is in contrast to most other existing methods which 
work either conditionally, or with an efficiency that is typically less than 1 for deterministic 
schemes |l9|. Moreover, the created quantum states have an average photon number much 
greater than one, thus are much brighter than those from other sources of non-Gaussian 
light such as single photon emitters. 

In the work of Schwendimann et ai, a master equation for exciton-polariton condensates 
was derived, where the "system" is considered to be the k = condensate, and the "reser- 



voir" are all the modes with k ^ [27|, |28[. The equations are essentially a variation of 
standard lasing equations, where there are a combination of one and two polariton loss and 
gain terms, which under suitable conditions gives rise to a condensate population. In Ref. 



29| the authors obtain the Fokker-Planck equations for the Wigner functions, and solve the 



equations numerically and analytically for the approximate case. The analytical form of the 
radial dependence of the Wigner function showed no trace of negativity, thus based on this 
alone it appears that exciton-polaritons cannot be used to generate coherent non-Gaussian 
light that would be useful for continuous variables. 

There is, however, an aspect which was not discussed in Ref. 29|, which generally 
applies to all lasing master equations. For simplicity let us first consider the standard lasing 
equations, with a hypothetical built-in non-linearity which will be clear that it is analogous 



to the polariton case. The Scully-Lamb lasing equation 



30| with such a non-linearity reads 



f^=-l[H....,p]-^C[a\p]-lc[a,p] 

H [p(aa^)^ + 3aa^ paa^ — Aa) paa)a + H.c] (1) 
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FIG. 1: The steady-state Wigner function for (a) a laser with no phase fixing K = (b) an injection 
locked laser with non-linear interaction U = 0.1 and K = 50. Parameters used are A = 3, B = 1, 
7 = 1. The negativities N = J cPa[W{a) — \W{a)\]/2 are marked on the plots. 



where 



U t t 



(2) 



is the non-hnear interaction term, and 



C[0, p] = pO^O + O^Op - 20 pO^ 



(3) 



is a Lindblad loss or gain term, with O an arbitrary operator. Here A is the gain coefficient, 
B is the self-saturation coefficient, 7 is the cavity loss rate, U is the non-linear Kerr inter- 



action (self- interaction for condensates). A similar system 



nonlinearity as we consider here, was analyzed in Refs. 
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'or a as opposed to the x^^^ 



The steady-state solution of these equations can be obtained by setting 



dt 



0. From 



the structure of the equation ([T]) it can be deduced that in the steady-state limit only the 
diagonal components are non-zero {s^. This is most easily seen by calculating the behavior 
of where Pk = ^nPn,n+k- For diagonal terms ^ = 0, so that Pq = 1 for all times. 
In contrast, the off-diagonal terms decay away, and in the steady-state limit tend to zero. 
Physically, this can be interpreted as the effect of phase diffusion, where any initial coherence 
becomes randomized in the long-time limit. Far above threshold, the density matrix of such 
a state is 

p=— d9\re'^){re'^\ = e-''' V (4) 



n=0 



where a = re*^ and \a) is a coherent state 



35-]. The state is diagonal in the density matrix. 



and the Wigner function has a "doughnut" shape with equal amplitude but a mixture of all 
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phases (see Figured^). The point is then that for the diagonal components of the density 
matrix, the non-hnear term H^^^ is completely ineffective, no matter how large U may be. 
This can be seen by looking at the matrix elements 

{n\[H,^,, p\\m) = ^ [n{n - 1) - m(m - 1)] Pnm, (5) 

which disappear for diagonal components. Since diagonal components only couple to other 
diagonal components in ([T]), the distribution is independent of U. 

How then is it possible to make use of the non-linearity? Clearly we must recover the off- 
diagonal terms in the density matrix, instead of the completely phase- diffused state (jlj). This 



can be achieved by employin g a s uitable phase fixing method such as injection locking 
feedback phase stabilization 



or 



371 . |38| . The effect of such phase fixing is to take a completely 
mixed state such as (jlj), which is immune to the interaction H^^^, and make it approach a 
pure state depending upon the strength of the phase fixing. For a pure state it is well known 
ihat the interaction Hi^t produces highly non-Gaussian states such as Schrodinger cat states 



4l|, and hence we may expect that the interaction will produce non-Gaussian features to 
the state of the light. 

For simplicity, in this paper we consider injection locking to be the phase fixing mecha- 
nism. This corresponds to adding a coherent term {sgI 

= ihKo (6a^e-*("-"")* - afe^e'^"""")*) (6) 

to the master equation ([T]), so that the first term is — ^[i/int -|- iJi^^^jP]- Here loq is the 
frequency of the (slave) laser, and we work in the rotating frame a = ase~^^°^, where as is 
the operator in the stationary frame, b is the annihilation operator for the injected laser, and 
CO is its frequency which has been factored out. We assume that bright, on-resonant (wq = u) 
coherent light is used for the injected master laser, so that we may replace b — Be~^^ , where 
i? > is a c-number. Assuming without loss of generality that the injection locking laser 
has phase ■?/' = 0, and converting ([T]) into a Fokker-Planck equation for the Wigner function. 
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we obtain 



dW 
Hi 



-K cosO- 
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(7) 



where K = KqB. The steady-state solution of such an equation is shown in Figure [T]d. We 
see that the phase fixing combined with an interaction U leads to a steady-state Wigner 
function with non-Gaussian features, signaled by the presence of a negative region. This 
negativity is very stable with respect to initial conditions and remains a feature of the 
distribution once the steady-state regime is reached. Of course, such a laser with a built in 
non-linearity will in practice have exceedingly small values of U, creating extremely small 
negativities in the Wigner function. We shall, however, see that for exciton-polaritons such 
non-linearities are within experimental reach. 

We first describe our proposed device configuration (Figure EJ^a)). A standard microcavity 
quantum well system is used, where exciton-polaritons are pumped by an excitation laser 
or electrical injection 42|. A secondary injection locking laser is introduced at normal 
incidence of variable intensity. Metal deposition on the surface of the microcavity allows for 
a trapping potential which confines the condensate area A. Metal deposition is known to 
create an effective potential which traps the exciton-polariton condensate |25]. The emerging 



photoluminescence is then used as the output light with non-Gaussian 



To describe this system, our starting point is the theory of Refs. 27 



properties. 

where a master 



equation for the polariton condensate is derived. Including the injection locking terms, the 
density matrix of the condensate evolves in time according to 



^ = -4[-f^int + -f^iock, p] - T2C[aa, p] - A2C[a^a\ p] 
dt n 

-(ri + 7o)£[a,p]-Ai/:[at,p]. (8) 

Here a is the annihilation operator for a A; = polariton, and Uq is its frequency (which 
may include a mean-field energy shift due to interactions). Fi 2 are one and two polariton 
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FIG. 2: (a) The schematic device configuration. A standard exciton-polariton microcavity with 
distributed Bragg reflectors (DBRs) is pumped to form a condensate in the quantum wells (QWs). 
Here we show the example of a AlAs/AlGaAs structure that could be used for the device, although 
other materials could be used. An injection locking master laser is injected perpendicularly in 
order to add a phase fixing term to the condensate. The results of this paper are concerned with 
light labeled as "output" emerging from the microcavity. (b)(c) Typical steady-state output of 
the injection locked polariton system. Parameters used are D2/D1 = 10, Gi/Di = 1, G2 = 0, 
U/Di = 10, K/Di = 10^ 



loss terms due to polariton scattering, Ai^2 are one and two polariton gain terms due to 
scattering, and 70 is a loss rate due to leakage through the microcavity mirrors. From ([8]) it 



is possible to derive a Fokker-Planck equation for the Wigner function 
techniques |39| we obtain 



29| . Using standard 



dW 



r or n 
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U fl d'^W d^W 1 d^W 
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where we have retained only derivatives up to second order in the Lindblad terms in (|8]), 
since this gives a very small correction to the overall dynamics. Here Gi = Ai — Fi, 



G2 = ^2 — r2 — 7o, -Di = Ai + Fi, D2 = A2 + r2 + 7o are the net gain and diffusion 
coefficients. The precise expressions for Gi and Di may be calculated using the expressions 



m 



27 



29| . While different materials give rise to different parameters, the integral expressions 
given in the above references possess certain properties which are common to all systems. As 
can be deduced from dimensional analysis, the order of the parameters are approximately 

At Ti Q Ao To (9) 

where Vq = Qe^as/ '^T^eA is the polariton-polariton interaction energy, A is the quantization 
area, a is a temperature length scale set according to Eq = ti^ /2m^^^a^ = ksT, where 
T is the temperature, and is the exciton Bohr radius. We can observe that the order 
0(Ai/A2) ~ 0(ri/r2) ~ 27r^a^/A is typically a small parameter such that Ai, Fi ^ A2, r2. 
However, Ai,ri cannot be neglected as without these condensation is unstable. Meanwhile 



the interaction can be estimated to be 40(] 



TT' 



where X is the excitonic Hopfield coefficient. From the expressions in Ref. [27|] it can be 

deduced that Gi < and Di 2 > 0, but Go can change sign which can be taken to be the 

2 g g 

criterion for condensation assuming Ai ^ (72 |27| . 

In Figure Mjo) and (c) we plot the typical Wigner distribution for the steady-state oper- 
ation of the injection locked polariton BEC. The distribution typically consists of a bright 
peak that is phase fixed due to the injection locking. Adjacent to the peak, there is a region 
of negativity, approximately in the direction of origin. The negativity appears as soon as 
the injection locking is switched on, and at first grows with both injection locking strength 
K and interaction strength U. For all parameters K > 0,U > we see stable negativities 
of the Wigner function, quantified by = / (Pa[W{a) — \W{a)\]/2. Figure |3^ shows the 
behavior of the amount of steady-state negativity with K. At large enough K, there is a 
saturation of the negativity, beyond which the total negativity no longer increases, and in 
fact decreases slowly beyond this point. Although the overall negativity decreases for very 
large K, the Wigner distribution tends to become more sharply defined for large K. For 
smaller values of K, the negative region tends to be shallower, but distributed over a larger 
area. In Figure [3]d we plot the negativities for D2/D1 = 10, which corresponds to very small 
(r = 0.16/im) traps or very high-Q cavities giving very long polariton lifetimes (I/70 = Ins). 
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FIG. 3: Dependence of the total negativity with injection locking strength K for various inter- 
actions U (in units of Di) for the polariton BEC. Steady-state values are taken for (a) estimated 
experimental parameters D2/D1 = 100, (b) very high Q cavity/small trap case D2/D1 = 10. Other 
parameters are Gi/Di = 1, G2 = 0. 

We see that the reduced phase diffusion improves the negativity of the Wigner function. We 
generally also observe that close to threshold gives the largest negativities. 

In conclusion, we have proposed a device for deterministically producing bright coherent 
non-Gaussian light from an exciton-polariton condensate. It was shown that steady-state 
generation of such light is possible using currently accessible experimental parameters when 
a suitable method of phase fixing is imposed on the polariton condensate. Although we 
have considered injection locking in this paper, other methods such as feedback control of 
the phase should accomplish the same task. The phase fixing overcomes the phase diffusion 
which tends to dephase the condensate and leads to a completely mixed state; with phase 
fixing, non-Gaussian states are created exploiting the natural nonlinearity of the polaritons. 
The completely deterministic nature of the production of these non-Gaussian states, as 
well as their brightness, is potentially useful in the context of continuous-variable quantum 
information processing, for which at least some resource states must possess a negative 
Wigner function. On-demand sources will save quantum memories and help scalability, while 
mean photon numbers above unity allow for making use of larger Hilbert spaces encoded 
into the infinite-dimensional space of an optical mode. 
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